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Abstract 

In this paper, several aspects of the turbulent relative particle dispersion are considered, 
many of which are motivated by laboratory experiment and numerical simulation results. 
These are, 

* spatial intermittency effects, 

(a) reduction of relative particle dispersion in 3D; 

(b) prevalence of power-law scaling of relative particle dispersion in 2D enstrophy cas- 
cade; 

* quasi-geostrophic aspects exhibiting an enhanced relative particle dispersion in the baro- 
clinic regime and a negative eddy- viscosity development for some insight into this aspect; 

* reduction of relative particle dispersion and prevalence of the ballistic regime due to 
compressibility effects. 
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1. Introduction 



Taylor [T] introduced the concept of diffusion by continuous movements and defined a 
diffusion coefficient D such that the mean square displacement {[s(t)} 2 ) is given by 

<[ S (t)] 2 > = 2^ (1) 

t being the elapsed time. Richardson [2], on the other hand, proposed that turbulent diffusion 
should be characterized by the distance between neighboring particles, because if the inter- 
particle distance is within the inertial range, one may expect to find universal super- diffusive 
behavior in the relative particle dispersion process, which may be interpreted in terms of an 
interparticle separation dependent turbulent diffusivityQ. From a purely empirical analysis 
of atmospheric dispersion data, Richardson [2] then showed that the turbulent diffusivity 
defined by the rate of increase of the mean square interparticle separation distance, 

Bs 5iwn (2) 

goes like 4/3 power of this distanced, 

D~([ J R(t)] 2 > 2/3 . (3) 

The most important application of two-particle dispersion is in understanding the motion of 
passive scalars (like pollutants) in atmospheres and oceans. 



1 Indeed, as Taylor '■'> mentioned, "Richardson was a very interesting and original character who seldom thought on the same 
lines as his contemporaries and often was not understood by them." 

2 On theoretical grounds, (3) can be justified only if the advecting velocity field is finite correlated in time. This follows by 
noting, from the initial-value problem 



f t R(t) = V(t) 
t = : R = Ro 



that 

2 

Assuming a stationary process, this gives 



1 -([R(t)-R ] 2 ) = J J (V(t')V(t"))dt'dt". 
oo 

t 

(\R(t) -Ro} 2 ) =2tj C{r)dr 



where C(r) is the autocorrelation, 

C(r) = (V(O)V(r)) . 

OO 

If C(t) decays in time fast enough so J C(r)dr is finite, we have a diffusive motion with turbulent diffusivity, 



oo 

v = t ^^ t (lR(t)-Ro] 2 )=fc(T) d T. 

o 

On the other hand, if C(r) decays slowly in time, we have a super-diffusive motion with 

([R(t) -Ro] 2 ) ~t 2 \ u>l/2 

and 

2v — l 

D ~ ([R(t) - Ro] 2 ) 2 " . 
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Obukhov |4j showed that Richardson's relation (3) can be derived via Kolmogorov's [5] 
(K41) theorjo for homogeneous isotropic 3D fully developed turbulence (FDT). When the 
interparticle separation is within the inertial range, Obukhov |3] gave 

D ~ e 1/3 ([i?(t)] 2 ) 2/3 (4) 

e being the mean energy dissipation rat^l, hence the name Richardson-Obukhov (RO) scaling. 
It may be mentioned that Richardson [2j indicated that D has the dimensions of e 1//3 and 
that T), like e, remains nearly constant. So, Richardson [2], interestingly provided the first 
experimental evidence for dissipative anomaly (i.e., e is almost independent of Reynolds 
number) and came very close to fully discovering the basic tenets of the Kolmogorov [5] theory 
15 years earlier. On the other hand, the Richardson formulation connects with the universal 
aspects of FDT actually stronger than the Kolmogorov formulation since the unphysical 
effects due to sweeping by large scales are precluded from the outset in the Richardson 
formulatiorj!. 

In this paper, several aspects of the turbulent relative particle dispersion are considered, 
many of which are motivated by laboratory experiment and numerical simulation results. 
These are, 

* spatial intermittency effects, 

(a) reduction of relative particle dispersion in 3D, corroborating the numerical simula- 
tion results (Boffetta and Sokolov [7]); 

(b) prevalence of power law scaling of relative particle dispersion in 2D enstrophy cas- 
cade, corroborating the difficulty in observing Lin [5] scaling law in laboratory 
experiments (Jullien [H]); 

* quasi-geostrophic aspects exhibiting some interesting departures from classical 2D re- 
sults and a negative eddy-viscosity development for some insight into these aspects; 

* particle clustering due to compressibility effects, corroborating the laboratory experi- 
ment and numerical simulation results (Cressman et al. [TO]) where compressibility is 
effectively produced on the free surface of a shallow fluid layer. 



2. 3D Relative Particle Dispersion 

Consider the temporal evolution of the separation between two particle trajectories in 
the inertial regime < C I in high R e FDT. The equation for the separation distance 
between two particles is 

jR(t) = 5v(R,t) (1) 



3 The Kolmogorov [5] theory stipulates that the probability distribution function of velocity fluctuations, in the inertial range, 
depends only on the mean energy dissipation rate e and the length scale i. Further, e is independent of the Reynolds number 

Re- 

4 There is some logical inconsistency in the expression in (4) - D and (^[R(t)]' 2 ^ are Lagrangian quantities while e is defined 
fully precisely in the Eulerian formulation. 

5 Indeed, as Monin and Yaglom [6] remarked, Richardson's formulation "indicates Richardson's belief in the existence of a 
universal physical law of sufficiently simple form." 



3 



from which, 



|[R(t)] 2 = 2R(t).5v(R,t) (2) 



Assuming the scaling relation, 

|5v(R,t)| ~ [R(t)] a , a<l (3) 

and averaging over many particles and many different initial separations, we obtain from 

(2), 

j t ([R(t)]) 2 -([R(t)} 2 f +a)/2 (4) 

from which, 

([R(t)] 2 ) ~ . (5) 

Observe that (5) shows a super- diffusive growth since a < 1 (for ordinary diffusion like the 
Brownian motion, ([R(t)} 2 } ~ t). For the K41 scaling given by a = 1/3, (5) leads to the RO 
scaling]! 

<[i?(t)] 2 )~t 3 . (6) 



3. 3D Relative Particle Dispersion: Effects of Spatial Intermittency 

The RO theory does not take into account the spatial intermittency in FDT that was 
revealed by laboratory experiments and numerical simulations. As is well known, spatial 
intermittency effects would cause systematic departures from the RO scaling law (6) which 
uses mean energy transfer rate. On the other hand, spatial intermittency effects are known 
to become more pronounced at small scales, so the turbulence activity gets concentrated in 
smaller and smaller regions of space and the active region (called the dissipative structures) 
becomes strongly convoluted like a fractal (Mandelbrot [H]). Indeed, the smaller the particle 
separation, the stronger the spatial intermittency effects become. The fractal aspects of FDT 
may be simulated in a first approximation by representing the dissipative structures via a 
homogeneous fractal with non- integer Hansdorff dimension D (Frisch et al. p2]). Using then 
the result, 

1 D -3 

« = o ^ o (') 



(5) becomes, 



([R(t)] 2 )~t 



3-3 



' 3- Dp 
J'- D 



Noting that D < 3, (8) shows that the effect of spatial intermittency is to cause reduction 
in the relative dispersion, in agreement with the numerical results of Boffetta and Sokolov 
[7]. Further, (8) also shows that the spatial intermittency effects, no matter how strong 



6 For a chaotic system (with positive Liapunov exponent A), we have 

R(t) ~ e xt . 



In FDT, there is a unique Sv (R, t) for each R, so there is now a continuum of \(R). Consequently, R(t) grows algebraically in 
FDT rather than exponentially. 
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(i.e. even in the limit Dq =>- 0), cannot change the super-diffusive nature of relative particle 
dispersion, i.e., 

([J2(f)] a > ~ f , 1/ > 1. (9) 

It is of interest to note that the relative dispersion of particle pairs in FDT was also 
considered by Roberts [13] using the Eulerian DIA method (and eliminating the spurious 
effect of sweeping motion of the small eddies by the energy-containing eddies). If the energy 
spectrum is given by 

E(k) ~ k-' n (10) 

Roberts [13] gave 

([R(t)} 2 ) (11) 
In the presence of spatial intermittency, we have from (7), 

n = I + I (3 - D ) . (12) 

Using (12), (11) gives 

([R(t)} 2 ) ~J 3+3 (^)] (13) 

which gives a relative particle dispersion faster than that given by the RO scaling law (6) 
contrary to the expected slow down of relative-particle dispersion by spatial intermittency. 



4. 2D Relative Particle Dispersion 

In 2D FDT, for pair separations larger than the energy injection scale but smaller than 
the integral scale (as in the inverse energy cascade), we have the scaling behavior, 

5v(R) ~ R 1/3 (14) 

which leads to the following growth relation, 

j t ([R(t)f)^([R(t)]Y 3 (15) 

and hence, the mean square pair separation grows as (Kowaleski and Peskin [14"]). 

([R(t)} 2 )~t\ (6) 

On the other hand, when the pair separation exceeds the integral scale, the particles become 
uncorrelated and Brownian diffusion sets in with the following growth relation, 

([R(t)} 2 )-t. (16) 

Okubo [15] analyzed oceanic experimental data and found the i? 4//3 law for the relative 
diffusivity along with the t 3 law for the mean square pair separation, in agreement with (15) 
and (6), respectively. The laboratory experiments of Jullien et al. [TB] also confirmed (6). 

For pair separations much smaller than the energy-injection scale, the velocity field would 
be smooth. The relative particle dispersion in the 2D enstrophy cascade was considered by 
Lin [8] who gave for the relative diffusivity, the following scaling result, 

j t ([R(t)f)~T^([R(t)f) (17) 
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and hence, for the mean square separation, the exponential growth behavior, 

([R(t)} 2 ) ~ e' 1/3 « (18) 

where r is the mean enstrophy dissipation rate. 

Balloon measurements in the atmosphere (Morel and Larcheveque [17] and Er-El and 
Peskin [18]) provided evidence comprising second-order structure functions and relative dif- 
fusivities which were approximately proportional to the square of the separation length and 
separation variances which grew approximately exponentially in time, hence supporting the 
existence of the enstrophy cascade and the Lin relative dispersion scaling result (18) in 2D. 
Babiano et al. [19] did numerical calculation of relative dispersion in 2D FDT and found 
that, if the initial pair separation is larger than the energy-injection scale the relative disper- 
sion follows the t 3 law (6) up to the most energetic scales. On the other hand, the relative 
particle dispersion was found to follow the e rl/3t law (18) only in a very short transient 
stage when the initial separation lies at the bottom of the enstrophy cascade. This was also 
confirmed by Kowaleski and Peskin [T3] via numerical calculations. Particle-path date from 
a float experiment (Ollitrault et al. [20]) showed that the relative diffusivity varies as R 2 
and -R 4 / 3 for distances smaller and larger, respectively, than a forcing scale of the order of 
the Rossby radius of deformation. Two particles initially separated between 40 km and 300 
km dispersed according to the t 3 law (6) (with relative diffusivity ~ i? 4//3 ) while those with 
smaller initial separation distances dispersed according to the e 7 " 1 3 * law (18) (with relative 
diffusivity ~ R 2 ). The scaling law (18) was confirmed at early times in a laboratory exper- 
iment by Jullien [9], but the range of pair separation scales over which this was observed 
was again very small, so this observation was rather difficult. This situation appears to be 
traceable to spatial intermittency effects in the enstropy dfi as shown in the following. 
It is of interest to note that even the numerical investigations of Elhmaidi et al. |30j and 
Zouari and Babiano [3T] showed that relative dispersion in 2D energy cascade is affected 
by coherent structures and indeed exhibits steeper than the t 3 law, symptomatic of spatial 
intermittency corrections. 



4.1 Inverse Energy Cascade 

The spatial intermittency effects in the inverse energy cascade may again be incorporated 
via the fractal aspects of active turbulence regions. Thus, using the result (Frisch et al. [12]). 

3 -A) , 1Q , 

a = — - — (19) 

(5) becomes, 

([R(t)] 2 ) ~ tl 3+ ^ 2 -M . (20) 

Noting that Dq < 2, (20) shows that the effect of spatial intermittency in the energy cascade 
is to make the relative particle dispersion go steeeper than the t 3 law, in agreement with the 



7 Direct numerical simulations of freely decaying 2D FDT (McWilliams 21, Benzi et al. [22], Brachet et al. [23], Kida [24], 
Ohkitani 1251 , Schneider and Farge |26l ) and forced-dissipative 2D FDT (Basdeveant et al. I27| . Legras et al. 1281 , Tsang et al. 
\ '29\ ) showed spatial intermittency in the enstrophy dissipation field caused by the presence of coherent structures. 
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numerical calculations of Elhmaidi et al. [30] and Zouari and Babiano [31 



4.2 Enstrophy Cascade 

On incorporating the spatial intermittency effects in the enstrophy cascade as per the 
homogeneous fractal model for the enstrophy dissipative structures, we have (Shivamoggi 

a = — — • ( 21 ) 

Using (21), (5) becomes 

([R(t)] 2 ) ~J 3+3 (^fe)]. (22) 

The effect of spatial intermittency is to make the relative particle dispersion go slower than 
the exponential growth law (18) which becomes operational in the space-filling limit Dq =>- 2. 
On the other hand, this also shows that, in the presence of spatial intermittency, however 
small, the exponential growth law (18) is replaced by the power law growth (22). This ap- 
pears to support the difficulty in observing the exponential growth regime in both laboratory 
experiments and numerical simulations. 



5. Quasi-geostrophic Relative Dispersion 



The dynamics of a 3D rapidly rotating fluid is characterized by the geostrophic balance 
between the Coriolis force and pressure gradient transverse to the axis of rotation. Quasi- 
geostrophic dynamics refers to the nonlinear dynamics governed by the first-order departure 
from this linear balance and is inherently 3D. The governing equation is the quasi-geostrophic 
potential vorticity equation (Charney [33J ) for an equivalent barotropic fluid in the /-plane 
(/ being the local Coriolis parameter). The term representing baroclinic effects in the flow 
in this equation introduces a characteristic length scale, namely the Rossby radius i?o = 
y/gH/ f, into the problem (H being the depth of the ocean taken to be uniform and g 
being the acceleration due to gravity). Consequently, this problem exhibits some interesting 
departures from the properties of classical 2D turbulence (Shivamoggi |32j). 



5.1 Energy Cascade 

Upon incorporating the spatial intermittency effects in the energy cascade as per the 
homogeneous fractal model for the energy dissipative structures, we have (Shivamoggi [32J) 

^ , In » R 

a = { (23a, b) 



3^Do o < J R ( 



0- 



8 Indeed, relative particle disperison in the inverse energy cascade was found (Babiano et al. |19| . Jullien et al. 1161 ). unlike 
Brownian motion, not to be a progressive process, but rather involving sequences of quiet periods and sudden bursts. 
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Using (23), (5) becomes 



<[«(*)] Vr£i r ln>>R0 (24a, b) 



The effect of spatial intermittency in the energy cascade is again to make the relative particle 
dispersion go steeper than the inertial-range scaling laws, 

(wm < 25a - b » 

Observe that the relative particle dispersion in the baroclinic regime (£ n 3> Rq) is greatly 
increased due to enhanced vortex stretching in this regime. 

Further insight into the unusual aspects quasi-geostrophic relative dispersion problem in 
the inverse energy cascade may be obtained by considering the eddy viscosity development 
for quasi-geostrophic turbulence. 



5.2 Eddy Viscosity for the Inverse Energy Cascade 

Kraichnan [31] proposed that the spontaneous development and net energy gain of large- 
scale structures in the inverse energy cascade from small-scale turbulence can be described 
by a negative eddy viscosity (which is either introduced phenomenalogically or derived 
via closure approximations)!^. Following Kraichnan [31] , one may treat the eddy viscosity as 
constant and calculate it by balancing the net eddy- viscous gain with the energy flux rate e 
into the explicit scales, 

2k 2 v T E{k)dk = -e (26) 

k c being the cut-off wave number so the explicit scales are given by k < k c . The negative sign 
on the right in (26) arises from the fact that energy is flowing toward smaller wavenumbers 
(rather than larger wavenumbers). 

Using the energy spectra (Shivamoggi [32]), 

p/.x / c k e^k~^ , kR » 1 

m ^\c k e^R-^k-^,kR «l (27a ' b) 



Cfc being a constant. 

Using (27), the energy balance relation (26) gives 



~ < -H-V/a (k c R f 3 k~^ , k c R « 1 (28a ' b) 



9 It should be mentioned that the concept of eddy viscosity is not on strong grounds because the basic underlying idea that 
scales of motion of given size are acted on by smaller scales as if the latter were an augmentation of the equilibrium thermal 
agitation is not totally valid, thanks to the lack of clear separation of the two scale sizes. 
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which may be rewritten in the Leslie-Quarini |35j universal form, 




(29a, b) 



(29a) corresponds to the result for classical 2D turbulence given by Kraichnan [31] while 
(29b) corresponds to the baroclinic regime - observe the explicit appearance, as one would 
expect, of the baroclinic parameter 

a = k c R (30) 

in the baroclinic regime^!. (29) shows that the Leslie-Quarini for the eddy viscosity has 
a certain robustness to it (as also indicated previously by Shivamoggi and Hussaini [37]). 
Observe further that the higher turbulent transport!"] in the baroclinic regime indicated by 
(28b) and (29b) is consistent with an enhanced relative particle dispersion indicated by (25a). 

For some insight into an actual physical mechanism underlying the negative eddy viscosity, 
see Appendix. 



5.3 Enstrophy Cascade 

Upon incorporating the spatial intermittency effets in the enstrophy cascade as per the 
homogeneous fractal model for the enstrophy dissipative structures we have (Shivamoggi 

^ , in > Ro 




a = { (31a, b) 

Using (31), (5) becomes 



<[*<«)r> ' T i >>ii0 < 32a - b > 

The effect of spatial intermittency in the enstrophy cascade is again to make the relative 
particle dispersion go slower than the inertial-range scaling laws, 

Observe the particle clumping indicated in the baroclinic regime (t n ^> R ) because the 
divorticity sheets favored in the enstrophy cascade baroclinic regime are more likely to occur 
near vortex nulls (Shivamoggi |38j) and particle clumping is further favored to occur near 



10 A similar result occurs for eddy viscosity in a compressible turbulence (Shivamoggi and Hussaini I37p showing the explicit 
appearance of the Zakharov-Sagdeev 1361 compressibility parameter. 
11 This may be seen by rewriting (28) in the following form 

UT ~ \_| Cfe -l £ l/3 i?0 4/3 a -2/3 _ a < L 
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vortex nulls (see Section 6). On the other hand, Charney [33] showed that, for £ n <C R , 
a potential enstrophy inertial range exists, in which the relative particle dispersion grows 
exponentially, as indicated by (33b). 



6. Compressibility Effects on Relative Particle Dispersion 

Intuitively, fluid compressibility is believed to lead to trapping of particles for long times 
and counteracting their tendency to drift away from each other - strong fluid compressibility 
would lead to particle clustering (Falkovich et al. [ID]). Laboratory experiments and nu- 
merical simulations in full-fledged 3D compressible FDT are not at hand yet. Shallow fluid 
layer flows provide an interesting alternative in this regard because the horizontal divergence 
of the free-surface flow on a shallow fluid layer is non-zero even though the fluid is incom- 
pressible. Consequently, the motion of passive tracer particles^ (used as surface markers 
in oceans) is not representative of 2D incompressible turbulence because such particles will 
respond only to the fluid flow on the free surface and not to the flow normal to the surface. 
Hence, they sample the horizontal components of the flow-velocity field and provide a con- 
venient framework to analyze the compressibility effects on relative particle dispersion even 
though the flow velocity is very small compared with the speed of sound (Sommerer and Ott 
[41]). Laboratory experiments and numerical simulations have therefore investigated relative 
particle dispersion on a free surface (Cressman et al. [10J). On the other hand, Cressman et 
al. [TO] expressed the necessity to have a theoretical framework to explain their results. We 
will now proceed to provide one such theoretical formulation. 

On assuming barotropic fluid and adiabatic flow processes, we obtain (Shivamoggi [42J) 



where the polytrope exponent 7 (1 < 7 < 00) may be treated as a compressibility parameter 
(the limit 7 =>■ 00 corresponds to the incompressible fluid limit). 
Using (34), (5) becomes 



(35) shows that the effect of compressibility is to make the relative particle dispersion exhibit 
a power law growth with a scaling exponent smaller than 3 (given by (6)), in agreement with 
the laboratory experiments and numerical simulations (Cressman et al. [TO])- 

On the other hand, on incorporating the spatial intermittency effects in compressible 
turbulence (as indicated by the numerical simulations of Lele et al. [43] Passot et al. [44]), 
and using a homogeneous fractal model for the kinetic energy dissipative structures, we have 
(Shivamoggi [45]). 



1 Particle tracers which have the same density as that of the carrier fluid and very small size can be approximated as point-like 
particles having the same velocity as that of the carrier fluid at the position of the particle. 




(35) 




(36) 
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Using (36), (5) becomes 



3-1/7 



(3-D ) 



(37) 



Noting that Dq < 3, (37) shows that the effect of spatial intermittency is to cause further 
reduction in the relative particle dispersion. On the other hand, on noting that the dissipative 
structures in a compressible turbulence are typically shock- wave like (D = 2), (37) becomes 



indicating that relative particle dispersion in intermittent compressible turbulence occurs 
in the ballistic regime. This appears to be physically plausible because, in the presence 
of shock waves, particle clustering^ renders the velocity increment become independent of 
the particle separation, and hence forcing the relative particle dispersion go ballistic. It is 
of interest to note, in comparison with (38), that laboratory experiments (Cressman et al. 
[TO] ) indicated a scaling exponent of 1.65 while numerical simulations (Cressman et al. [TO] ) 
indicated a scaling exponent of 1.80. Further, (37) also shows that compressibility effects, no 
matter how strong, cannot change the super- diffusive nature of relative particle dispersion. 
Indeed, in the infinite compressibility limit 7 1, (37) yields 



pertaining again to the ballistic regime (38)! 
7. Discussion 

In this paper, several aspects of the turbulent relative particle dispersion are considered, 
many of which are motivated by laboratory experiment and numerical simulation results. 
These are, 

* spatial intermittency effects, 

(a) reduction of relative particle dispersion in 3D, corroborating the numerical simula- 
tion results (Boffetta and Sokolov [7]); 

(b) prevalence of power-law scaling of relative particle dispersion in 2D enstrophy cas- 
cade (no matter how weak spatial intermittency effects are), corroborating the 
difficulty in observing Lin [8] scaling law in laboratory experiments (Jullien [9]) 
and numerical simulations (Babiano et al. [19] and Kowaleski and Peskin [Ti]); 

* quasi-geostrophic aspects exhibiting an enhanced relative particle dispersion in the baro- 
clinic regime and some insight into this aspect has been explored via a negative eddy 
viscosity development; 



13 Thanks to a comprossiblc-fiuid flow like flow situation prevailing on the free-surface of shallow fluid layers, local convergence 
and divergence regions of particle density are observed there (Elhmaidi et al. 1301 and Cielsik et al. 46 ). In convergence 
regions where particles tend to clump there is downwelling and vice versa. Further, downwcllings are found to occur near 
strain-dominated (thin elongated) regions while upwellings are found to occur near rotation-dominated (patch like) regions. 



([ J R(t)] 2 )~t 2 ,V 7 





(39) 
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* reduction of relative particle dispersion and prevalence of the ballistic regime due to 
compressibility effects, corroborating the laboratory experiment and numerical simula- 
tion results (Cressman et al. [10]) where compressibility is effectively produced on the 
free surface of a shallow fluid layer. 

It should be mentioned, however, that the 3D RO scaling result (3) has received little 
experimental support due to the difficulty of performing Lagrangian measurements over a 
broad enough range of time and with sufficient accuracy. Even in recent laboratory experi- 
ments (Ott and Mann [17], Sawford [H], Bourgoin et al. [49J, Salazar and Collins [50]) with 
high-speed photography to track particles, and in numerical simulations (Yeung [51]) with 
the highest possible resolution possible for homogeneous isotropic turbulence, it is known to 
be hard to obtain an extended range with the RO scaling. The difficulty appears to be due 
to, 

* contamination of the inertial range by dissipative effects at the ultraviolet end and the 
external forcing effects at the infrared end of the spectrum caused by inadequate scale 
separation; 

* persistent memory of initial separation. 
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Appendix. Physical Phenomenological Derivation of Negative Eddy 
Viscosity 

A sound phenomenological demonstration of negative eddy viscosity via a simple model 
that captures the essential physics is apparently not at hand despite Kraichnan's [31] at- 
tempts motivated by complicated analytic developments. Analytical approaches used by 
Sivashinsky et al. [22], [S3] and Gama et al. [51] are restricted to flows possessing special 
symmetries. We propose to contribute toward this quest and give a simple phenomenological 
derivation of negative eddy viscosity for 2D energy cascade. 

We borrow from a basic idea for a physical mechanism underlying the negative eddy 
viscosity propsed by Kraichnan [33] - the large scales strain the small scales while a secondary 
flow associated with small scales grows. Further, the deterministic large-scale flow and a 
random small-scale vorticity field are characterized by disparate scales (Tur et al.|55j and 
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Fidutenko [56]). Consider therefore a random homogeneous small-scale vorticity field uj 
superposed on a stationary large-scale flow with stream function with velocity given by 

V = i 2 xVf (A-l) 

Linearizing about the large-scale flow (in analogy with a phenomenological development for 
MHD turbulence sketched by Pouquet [57J) we obtain, for the small-scale flow, with velocity 
v, 

^ = uV* (A-2) 
dt 

where, 

u> = V x v = oui z . 

Assuming a normal-mode type evolution in time for the small-scale flow, 

v ~ e at (A-3) 
so 1/cr may be interpreted as a coherence time of the small-scale flow, equation (A-2) becomes 

v = -W. (A-4) 
a 

On the other hand, the large-scale vorticity Q evolves according to 

^ (O) + V • (vw) = 0. (A-5) 



Assuming the large-scale flow to comply with the Beltrami condition 14 !. 

(0> = L 2 <*) (A-6) 
L being a constant (of dimension length), and using equation (A-4), equation (A-5) becomes, 

^W = 4W V2 * (A- 7) 

which implies that the nonlinear growth rate of the flow has the form of negative eddy 
viscosity given by, 

L 2 

v T (u 2 ) . (A-8) 



a 



Noting the spectral behaviors, 

(A-8) becomes 
in agreement with (28a)! 

14 A similar assumption was also made in the negative- viscosity development of Sivashinsky et al. |52| . |53| , 



v T ~ e 1/3 k- 4 ? 3 (A-10) 
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On the other hand, on recognizing that turbulent transport in the inverse cascade is ac- 
tually a competition between flow advection and vortex coalescence, (A-8) may be rewritten 

as 

v T ~^{{u 2 )-L 2 {u 2 )). (A-ll) 



A phenomenological derivation of the first term in (A-ll) may be given as follows. Con- 
sider a flow velocity v = (u, v) with a displacement vector (. = (£\, ii). If cr -1 is a typical time 
scale, we have on complying with the continuity equation (a la Batchelor [58] and Kraichnan 
PI) 

u ~ dli/dt ~ a£i 
v ~ dtijdt ~ — cr^2 



(A-12a, b) 



while the vorticity is given by 



v 11 /II 
w „___„_ ti £ a (_ + j. (A-13) 



The eddy viscosity (or turbulent diffusivity) is given by 

1 d 



v T ~--(e)~ a ((ll)-(lt)). (A-14) 



Using (A- 12) and (A-13), (A-14) becomes 



a 

in agreement with (A-ll). Here, 



v t ~-({u 2 )-L 2 {uj 2 )) (A-15) 



Observe that, according to (A-ll) or (A-15), the negative eddy-viscosity effect disappears 
if the small-scale flow decorrelates rapidly in time (i.e., a =^ oo), as Kraichnan |34j pointed 
out. 
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